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Abstract. We study a sample of 44 low-luminosity radio-loud AGN, which represent a range of nuclear radio-power 
spanning 5 orders of magnitude, to unveil the accretion mechanism in these galaxies. We estimate the accretion 
rate of gas associated with their hot coronae by analyzing archival Chandra data, to derive the deprojected density 
and temperature profiles in a spherical approximation. Measuring the jet power from the nuclear radio-luminosity, 
we find that the accretion power correlates hnearly with the jet power, with an efficiency of conversion from rest 
mass into jet power of ~ 0.012. These results strengthen and extend the validity of the results obtained by Allen 
and collaborators for 9 radio galaxies, indicating that hot gas accretion is the dominant process in FR I radio 
galaxies across their full range of radio-luminosity. 

We find that the different levels of nuclear activity are driven by global differences in the structure of the galactic 
hot coronae. A linear relation links the jet power with the host X-ray surface brightness. This implies that a 
substantial change in the jet power must be accompanied by a global change in its ISM properties, driven for 
example by a major merger. This correlation provides a simple widely applicable method to estimate the jet-power 
of a given object by observing the intensity of its host X-ray emission. 

To maintain the mass flow in the jet, the fraction of gas that crosses the Bondi radius reaching the accretion disk 
must be > 0.002. This implies that the radiative efficiency of the disk must be t; < 0.005, an indication that 
accretion in these objects occurs not only at a lower rate, but also at lower efficiency than in standard accretion 
disks. 
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1. Introduction. 

A fundamental problem in studying the physics of active 
galactic nuclei (AGN) identifying the mechanism that con- 
trols the level of activity. More generally, we need to de- 
termine why nuclear activity occurs in active galaxies but 
not in quiescent galaxies, and whether there is a disconti- 
nuity in the physical properties of these two galaxy classes. 
Addressing these points is also important in light of the 
prominent role that AGN feedback plays in the process of 
galactic evolution. Furthermore, AGN activity is related 
to the growth, via accretion, of supermassive black holes 
(SMBH); the clear evidence for a co-evolution of SMBH 
and host galaxies, represented by the connections between 
galaxy physical properties and SMBH mass, implies that 
the accretion process in AGN has a powerful influence on 
galaxy evolution. 

Several comparisons between the properties of active 
and quiescent galaxies have been performed with in at- 



tempts to isolate the quantities that correlate most sig- 
nificantly with the level of nuclear activity. We limit 
ourselves to radio-loud g alaxies in this paper. S ince the 
research of for example lAuriemma et al.l ( 19771) it has 
been known that the probability of an object to have 
a radio luminosity above a certain threshold is propor- 
tional to its optical luminosity. This was confirmed by 
studies of the bivariate radio/optical lu r ninosity function 
of ea rly-type galaxies ([Best et al ] 120051: iMauch fc Sadleil 
20071) . Objects of a given optical luminosity can be as- 



sociated with radio sources spanning a large range of ra- 
dio power. Several authors explored the role of environ- 
ment in setting the level of radio-emission. The fraction 
of radio-loud AGN is largely independent of the local 
galaxy density, but a depende nce is found on the larger- 
scale environment ( Best \2004 . A statistical trend links 
the amount of dust measured from the absorption fea- 
tures, seen in their o p tical images, and radiq-pow- er (e.g. 
de Ruiter et"al]|2002l) . IValentiin fc Biileveldl (|l983l ) found 
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a positive correlation between the X-ray luminosity of the 
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hosting cluster and the r adio- luminosi ty of the dominant 
cluster galaxy. Similarly. iBurnsI ( 199Clf ) showed that cen- 
tral, massive galaxies in clusters or groups of galaxies with 
a cooling flow are more likely to be radio-loud AGN than 
galaxies of similar mass that are not centered on a cooling 
flow, suggesting that concentrated cooling atmospheres 
can stimulate radio-loud AGN activity. These studies pro- 
vide statistically important results, but it is still impos- 
sible to ascertain from any of these measurements, on an 
object-by-object basis, whether a galaxy hosts an active 
nucleus. This implies that we cannot yet associate an ob- 
servable quantity and, consequently, physical process with 
the level of nuclear activity. 



2. Sample selection and data preparation 

In Capetti fc Balmaverdd ( 20051 ). we focused on a sam- 



Su bstantial progress was reported by lAllen et al 
( 2006f ). They studied 9 nearby, low-power radio galaxies, 
selecting objects with well-defined cavities in X-ray emit- 
ting gas, which were cospatial with the radio-lobes: these 
cavities were most likely inflated by radio-jets. They were 
able to estimate the jet, kinetic power required to inflate 
the cavities. On the other hand, it was possible, using the 
same X-ray images, to estimate the amount of hot gas 
available for accretion onto the SMBH. This analysis re- 
vealed a significant correlation between the accretion rate 
and jet power, which had a dispersion of only 0.16 dex. 
This represents not only an accurate and powerful method 
to predict the AGN energy output, but implies that spher- 
ical hot-gas, accretion, associated with the hot corona, is 
the main energy source powering the active nucleus. 



We extend the analysis performed by en et al 



to two samples of galaxies that host radio-loud 
nuclei of relatively low-power, which cover almost five 
orders of magnitude in radio-power. The comparison of 
accretion properties across the widest possible range of 
luminosity is clearly the next crucial step in unveiling 
the connections (and diversities) between galaxies show- 
ing different levels of nuclear activity. We examine data 
for the sample of 29 early - type "core" galaxies, selected 
bv iBalmaverde fc Capetti ( 2006f ). and for the sample of 
33 FR I radio galaxi es, which were ex t racted from the 3C 
sample as defined in IChiaberge et al. I (Il999l) : in the com- 
bined set of galaxies 44 have observational data in the 
Chandra archive (more details o f the sample select ion are 
given in Sect. [2]). With respect to Allen et al. ( 2006f ) we es- 
timated the jet power using a relation betwee n radio-core 
luminosity and jet kinetic power derived by iHeinz et al 



This relation provides a straightforward, robust 
measure of the jet power (with an rms error of 0.4 dex), 
based only on the core power, and enables the number 
of objects for which accretion and jet power can be com- 
pared, to be increased. Furthermore, the range of radio- 
power that can be studied covers almost five orders of 
magnitude; this corresponds to the full range of radio- 
luminosity of FR I, spanning from the faintest level of 
detectable activity in radio galaxies up to include galaxies 
whose radio-emission exceeds the luminosity marking the 
transition between FR I and FR II objects. 



pie of luminous, nearby early-type galaxies to study the 
link between their host properties and the active nucleus. 
We used archival HST observations to analyze their sur- 
face brightness profiles and to separate these early-type 
galaxies into core and power-law galaxies on the basis 
of the slope of their nuclear brightness profiles, follow- 
ing the Nuker scheme ( Lauer et al.lll995l ). We obtained a 
sub sample of 29 "core" galaxies (hereafter CoreG), that is 
galaxies for which the surface-brightness profiles were fit- 
ted with a double powe r -law o f innermost slope 7 < 0.3. In 
Balmaverde fc Capettil (|2006l ) we demonstrated that the 
CoreG hosted a radio-loud nucleus and that they could 
be considered to be miniature radio galaxies. In fact, they 
are drawn from the same population of early-type galax- 
ies as the FR I hosts, according to of the surface bright- 
ness profiles, and distributions of optical luminosities and 
black-hole masses, they differ from "classical" FR I ra- 
dio galaxies only in terms of their lower level of nuclear 
activity. 

We combine these CoreG with a sample of low- 
luminosity radio galaxies of morpholog y FR I extrac ted 
from the 3C cata logue of radio-sources _|Bennetll[l963 ) m 
a way defined bv IChiaberge et aL ( 19991 ). Three galaxies 
(namely 3C 270, 3C 272.1, and 3C 274) are in common 
between the two sample^ All 9 galaxies (wi t h the excep- 
tion of NGC 507) considered bv I Allen etall (|2006l ) are in 
common with our combined FR I/CoreG sample. 

We searched for Chandra observations available in the 
public archive up to August 2007 and found data for 22 of 
the 33 objects belonging to the FR I sample: w i th res pect 



to the sample considered bv IBalmaverde et"al ( 20061 ). 



we 

were able to add three sources observed, 3C 264, 3C 293, 
and 3C 442. W e found data for 24 of the 29 CoreG: 
with respect to ( Balmaverde fc Capettil 120061 newly ob- 
tained observations were found for UGC 0968, UGC 7797, 
UGC 8745. All together, we consider 44 galaxies that have 
a range of nuclear radio-power covering 5 orders of mag- 
nitude. The information on the Chandra observations and 
basic data for these galaxies are presented in Table [TJ 

We reduced all data using the Chandra Data analy- 
sis CIAO v3.4, with the CALDB version 3.3.0. We repro- 
cessed the data from level 1 to level 2 in the standard way, 
that is using the CIAO tool acis_process_events, removing 
bad pixels and applying pixel and PHA randomization. 
When possible, we improved the rejection of cosmic events 
for VFAINT mode observations. We rejected time inter- 
vals corresponding to high background levels and we fil- 
tered by selecting a standard set of ASCA grades (02346). 

We modeled the spectra obtained using XSPEC ver- 
sion 12.3.0. We used only events filtered between 0.3 and 8 
keV for our analysis and rebinned the spectrum to a min- 
imum of 20 counts per bin to ensure that a statistics 
could be applied. 



^ For simplicity, we include these sources in Table \T\ in the 
FR I sample. 
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Name 


Chandra 


summary 








Basic 


Data 






Obs. Id 


Exp. time 


D 


Log V 






Mk 


Log {Mbh/Mq) 


3C 028 


3233 


50.38 


666.7 


<38.86 


— 


-25.71 


8.99 


3C 031 


2147 


44.98 


67.15 


39 


40 


278 


-25.65 


8.70 


3C 066B 


828 


45.17 


83.68 


39 


90 


— 


-26.42 


9.32 


3C 075 


4181 


21.78 


91.0 


39 


30 


301 


-25.71 


8.84 


3C 078 


3128 


5.23 


112.1 


40 


88 


271 


-26.18 


8.61 


3C 083.1 


3237 


95.14 


98.65 


39 


10 


— 


-26.84 


9.51 


3C 084 


3404 


5.86 


69.31 


42 


10 


259 


-26.08 


8.58 


3C 189 


858 


8.26 


165.8 


40 


54 


— 


-26.22 


9.23 


3C 264 


4916 


38.33 


85.46 


39 


91 


271 


-25.17 


8.67 


3C 270 


834 


35.18 


29.69 


39 


22 


309 


-25.10 


8.72'' 


3C 272.1* 


803 


28.85 


14.11 


38 


57 


282 


-24.53 


9.00'' 


3C 274* 


1808 


14.17 


17.38 


39 


90 


333 


-25.39 


9.53'' 


3C 293 


5712 


3.10 


174.1 


40 


29 


185 


-25.36 


8.14 


3C 296 


3968 


50.08 


96.99 


39 


62 


299 


-26.17 


8.80 


3C 317 


890 


37.23 


134.3 


40 


64 


216 


-26.09 


8.26 


3C 338* 


497 


19.72 


118.7 


39 


96 


310 


-26.50 


8.89 


3C 346 


3129 


46.69 


569.3 


41 


74 


— 


-26.03 


9.14 


3C 348 


1625 


15.00 


544.9 


40 


36 


— 


-26.13 


9.19 


3C 438 


3967 


47.9 


907.7 


41 


14 


— 


-26.47 


9.34 


3C 442 


6392 


33.13 


103.1 


38 


12 


— 


-24.39 


8.40 


3C 449 


4057 


29.56 


67.46 


39 


06 


253 


-24.96 


8.54 


3C 465 


4816 


50.16 


118.1 


40 


37 


356 


-26.64 


9.14 


UGC 0968 


6778 


15.14 


32.39 


36 


94 


253 


-25.39 


8.54 


UGC 5902 


1587 


31.9 


12.68 


35 


83 


207 


-24.25 


8.00'' 


UGC 6297 


2073 


39.0 


13.63 


36 


46 


224 


-23.68 


8.33 


UGC 7203 


3995 


5.13 


31.89 


37 


44 


184 


-24.08 


7.98 


UGC 7386 


398 


1.43 


10.72 


38 


38 


238 


-22.97 


8.43 


UGC 7629* 


321 


40.1 


12.53 


37 


73 


291 


-25.09 


8.78 


UGC 7760* 


2072 


55.14 


5.227 


37 


30 


253 


-21.86 


8.54 


UGC 7797 


6785 


15.19 


29.77 


38 


05 


224 


-24.61 


8.33 


UGC 7878* 


323 


53.05 


14.83 


36 


90 


203 


-24.43 


8.16 


UGC 7898 


785 


37.35 


16.41 


37 


46 


336 


-25.34 


9.30'' 


UGC 8745 


6787 


15.18 


27.95 


37 


81 


232 


-25.07 


8.39 








25.28 


37 


31 


9'H8 
^oo 




o.'±o 


UGC 9723 


2879 


34.18 


13.79 


36 


92 


159 


-23.82 


7.73 


NGC 1316 


2022 


30.23 


20.64 


37 


82 


228 


-25.99 


8.36 


NGC 1399 


319 


56.66 


16.23 


37 


20 


344 


-24.75 


9.07 


NGC 3557 


3217 


371.99 


38.09 


37 


94 


273 


-25.70 


8.67 


NGC 4696* 


1560 


85.84 


36.89 


38 


65 


254 


-25.69 


8.55 


NGC 5128 


463 


19.6 


5.2 


39 


05 


120 


-24.64 


8.38'' 


NGC 5419 


4999 


15 


53.92 


38 


42 


332 


-26.14 


9.02 


IC 1459 


2196 


60.17 


19.97 


39 


38 


306 


-24.70 


9.18'' 


IC 4296 


3394 


25.4 


48.17 


39 


46 


336 


-25.69 


9.04 


NGC 507* 


2882 


44.21 


65.02 


36 


45 


310 


8.688 


8.89 



Basic data and Chandra observations log for the radio galaxies of the 3C/FR I and CoreG sample (with the addition of 
NGC 507). Objects from the Allen's sample (marked with *). Column description: (1) name, (2) Chandra observational 
identification number, (3) exposure time [ks], (4) distance [Mpc], (5) nuclear radio luminosity (5GHz) [erg s~^], (6) stellar 
velocity dispersion [km s~^] from the HyperLeda database, (7) total K band galaxy's absolute magnitude from 2MASS (except 
for 7 sources, see Sec. 13. 3p . (8) logarithm of black-hole mass in solar unity (see text for details). 



We adopt Ho=75 km s"^ Mpc~^ and qo=0.5. 



3. Calculation of the Bond! accretion rates 

Under the assumption of spherical symmetry and negligi- 
ble angular momentum, the accretion rate onto a SMBH 



can be estimated by the iBondil (|l952l ) approximation as: 
Mb ^ n Xcs pb tb^ 

where A takes the value of 0.25 for an adiabatic index 5/3, 
Cg is the sound speed in the medium (c^ — -y/F kT / p, nip), 
tb = '^'^^^fH ^jj^g Bondi accretion radius, and ps is the 
gas density at that radius, (p = l.lSneJTip), so that 
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Fig. 1. Example of the procedure followed to determine the Bondi power for NGC 7898: in the left 
panel, we show the spectrum extracted in four annuli, centered on the nucleus, that was modeled with a 
PROJECT*PHA(MEKAL+POWERLAW) law; in the right panel, we show the corresponding deprojected tempera- 
ture profile of the annuli. 



Mb = 



A.52ttXG\ 



The calculation of the Bondi accretion rate then re- 
quires the measurement of three quantities: the electron 
gas density at the Bondi radius, the gas temperature at 
the same radius, and the mass of the central black hole. 

Unfortunately not even the Chandra telescope, which 
provides presently the highest angular resolution in X- 
rays, is able to resolve the accretion radius in most galax- 
ies. Different strategies have been devised to estimate the 
gas and t emperature d e nsity at the accretion radius. For 
example, iGliozzi et al.l (j2003l ) fitted the electron density 
profile with a model of gas in i sothermal, hydrostatic e qui- 
librium in spherical symmetry. iHardcastle et al.l (|2007l ) de- 
rived the Bondi power for radio galaxies from the 3CRR 
catalogue, adopting density and temperature values com- 
parable to that typically measured (p = 5 10^ mT^ and 
kT = 0.7 keV) in nearby FR I. However, the results of 
Allen et al. ( 20061 ) indicated strong departures from equi- 



librium profiles and large variations in the central gas 
densities, suggesting that the accretion rates derived with 
these assumptions can be affected by large errors. 

In this work, we determine t he Bondi ac c retion rate 
following the same technique as en et al.l (|2006f ) and 
iBirzan et all ([iool . which is described in detail below. 
Briefly, we deproject the temperature and electron density 
profile, assuming spherical symmetry. We then extrapo- 
late the measured density to the accretion radius and de- 
scribe the density profile with a power-law, assuming that 



the temperature is constant inside the innermost measure- 
ment radius. Given the temperature, the electron density 
at the accretion radius, and the black-hole mass, we cal- 
culate the Bondi accretion. 



3.1. The IGM temperature profile 

To characterize the thermal and emissivity properties of 
the intergalactic hot gas, we extract the spectrum in four 
annuli (see Fig. [U left panel) centered on the peak of X-ray 
emission. Since a X-ray point source associated with the 
AGN is present in most galaxies, we extract the spectrum 
of the IG medium avoiding the innermost region within a 
radius of I'.'S. The outer radius is set to where the diffuse 
emission reaches the background level. 

To measure the dependence of temperature on radius, 
we deproject the spectra by assuming spherical symmetry 
and using the PROJECT model in XSPEC. This model 
projects prolate ellipsoidal shells onto elliptical annuli, cal- 
culating the projection matrix of each of the overlying 
shells onto the area of each annulus and computing the 
total projected emission on that annulus. The fitting oc- 
curs simultaneously for all annuli and, after completion, 
the model provides the best-fitting temperature and abun- 
dance values in each shell (Fig. [1] right panel). We adopt 
a single-temperature plasma model (MEKAL) and fore- 
ground absorption (PHABS), taking into account projec- 
tion effects, that is the model PROJECT*PHA(MEKAL). 
The hydrogen column density was tied between annuli and 
allowed to vary. With this fitting procedure we determine 
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UGC 7898 



1.0000 F 




Radius (pixel) 




0.00 



Radius [Kpc] 



Fig. 2. Example of the procedure followed in this work to determine the Bondi power for NGC 7898: in the left panel, 
we present the X-ray surface brightness profile of the galaxy; in the right panel, we show the derived deprojected 
electron density profile and the power-law extrapolation used to estimate its value at the Bondi radius (marked with 
a diamond). 



the deprojected profile of temperature. We then assume 
that the temperature measured inside the innermost an- 
nulus (usually bound between 1'.'5 and 5" in radius) is a 
good estimate of the temperature at the Bondi radius (see 
Table [2]). 

3.2. The electron density profile 

We determine the X-ray surface brightness profile (SBP) 
extracting, the counts in a series of circular annuli cen- 
tered on the central AGN (see Fig. [21 left panel, for an 
example; all SBP are shown in the Appendix A (Fig. I A. l]) . 
When present, we mask the X-ray jet and other unrelated 
sources; we then subtract the background measured far 
from the source of interest. To avoid contamination from 
non-thermal emission due to the central AGN, we exclude 
the inner 1.5 arcsec region. When a point source is present, 
we use the PSF libraries to create an image of the PSF for 
the specific observation. We normalize it to the source flux 
in a circle of 1 arcsec centered on the peak of the emission 
and subtract the wings of the PSF brightness profile from 
the profile of the source. 

As a general rule, we extract the X-ray counts from 4 
annuli in the first 5 arcsec, then increase the size of the 
extraction region outward. For weaker sources, we use a 
coarser sampling of the brightness profile to maintain suf- 
ficient signal-to-noise ratio. The next step is to deproject 
the observed SBP, to derive the number of counts emit- 
ted per unit volume as a function of radius. Assuming 
spherical symmetry, the counts contribution provided by 



each spherical shell to the inner ones can be determined by 
pure geome t ric con siderations, following the calculation by 



Kriss et al. 



Briefiy, we express the observed SBP 
as a matrix product between the deprojected count rate in 
the spherical shell and a weight matrix, defined as the por- 
tion of volumes of each overlying shell seen in projection 
on each annulus. 

Our aim is to compare the observed counts per unit 
volume, with the theoretical value expected for a ther- 
mal plasma of a specific electron density and metal abun- 
dances. The XSPEC fitting model provides us with the 
value of the thermal normalization and the model pre- 
dicted count rates (deprojected) in each of the annuli. 
Comparing these values with the counts per unit volume 
obtained by deprojecting the SBP, we adjust the thermal 
normalization to recover the counts emitted in each of the 
circular annuli used for the spatial analysis. The density 
was then calculated from the normalization of the ther- 
mal component, assuming Ue = 1.2 njy (for a fully-ionized 
gas) . The definition of the thermal normalization is 



K 



10 



-14 



47ri?^(l + z)2) 



where z is the redshift of the source. Da is the angular- 
diameter distance at that redshift, Ug and Tip are the elec- 
tron and proton number-density, respectively, and the in- 
tegration is performed over the volume of the projected 
shells. Assuming a uniform density distribution along the 
shell, we invert this equation to solve for the electron den- 
sity. Iterating this procedure for each spherical shell, we 
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obtain the electron density profile (see Fig. [21 right panel, 
for an example; all profiles are shown in the Appendix A, 
Fig.ESl). 

To estimate the gas density at the accretion radii, we 
perform a fit on the innermost points of the de nsity profile 
in a p ower-law form, Ueir) ~ r~", following lAllen et al 



3C83. 



( 2006f) and as explained in more detail in Sect. 13.51 



3.3. The black-hole mass 

The mass of the central black hole can be determined in 
various ways: when available, we prefer to use the BH mass 
estimates obtained fro m direct gas or s t ellar dynamical 
measurements listed in Marconi fc Hunt ( 20031 ). 

Otherwise, we derive the SMBH masses using 
the stellar velocity dispersion (from the HyperLeda 
database database). We apply the relation derived by 
Tremaine etal\ (|2002l ): 

log (MBff/Mo)=8.13 + 4.02 log (cr / 200 km s^i) 
which has an intrinsic dispersion of about 0.23 dex. 
Alternatively, for the 7 objects for which velocity disper- 
sion measurement is unavailable, we estimate the SMBH 
mass using the K-band absolute total magnitude (from 
2MASS data archivj^ by using the relation derived by 
Marconi fc mi^ (|20Cl3l) : 

log (Mbh/Mq) = 8.21 + 1.13 (log Lk -10.9). 
The observed scatter in this case is ~ 0.5 dex. 



3.4. Error budget 

The uncertainties in the results of the deprojection were 
estimated using a Monte Carlo simulations. We perturbed 
randomly each point of the SBP with a Gaussian distri- 
bution of the perturbations, which had an amplitude de- 
termined by the error bar. We obtained 1000 different re- 
alizations of the SBP. Each of them was deprojected, and 
we evaluated, in each shell, the mean value of the counts 
rate and the scatter. Then we propagate the errors in the 
calculation of the electron density, taking into account the 
errors in the counts rate and in the temperature measure- 
ment at the accretion radius. 

Another source of uncertainty is associated with the 
extrapolation of the electron density to r^. This is related 
to the uncertainties in the parameters of the power-law 
fit describing the density behavior but also to the choice 
of the range of radii to be included. We evaluated the 
differences in Ug, repeating the extrapolation using from 
3 to 6 points of the electron density profile. The adopted 
value for Ug is the average value of these measurements and 
its uncertainty spans the full range given by the overlap 
of each individual error bar. This method is graphically 
explained in Fig. [3l 



2 For 3C 066B, 3C 075, 3C 083.1, 3C 338, 3C 449, and 
3C 465 the K ban d luminosity is converted from H-band 
l|Donzelli et al.ll2007^ ■ assuming H - K = 0.21 and for 3C 442 w e 
convert the V-band luminosity from lSmith fc HeckmanI l|l989l) . 
assuming V - K — 3.30 



KOOO 



^ 0.010 



1 .0 

Radius [Kpc] 



10.0 



Fig. 3. Example of the procedure followed to determine 
the density at the Bondi radius in the case of 3C 83.1. 
The dotted (dashed-dotted) line represents the best-fit 
linear relation to the density profile obtained by using 
the innermost 3 (6) density measurements, while the dia- 
mond (star) is its extrapolation to the Bondi radius. The 
adopted value for Ug is the filled circle at r^. 



Concerning uncertainties in the black-hole mass deter- 
mination, we assume conservatively an error of a factor 3 
(0.5 dex), regardless of the method used. 

The final error in Pb is estimated by propagating the 
errors of all relevant quantities. We note that since the 
Bondi accretion power is tied to the electron density, re- 
lated in turn to radius via the value of a, the effective 
dependence of Pb on the black-hole mass is not simply 
from M^^, but from M^^. Given that the mean a value 
estimated (see Tabled]) is 1, the propagated uncertainty 
is essentially linear with the black-hole mass error. 

3.5. Accretion rate estimates and data quality 

The methods presented in the previous sections were ap- 
plied to every source of our sample to estimate the accre- 
tion power. However, the quality of the data is clearly not 
always sufficient to obtain a measure of Pb at a level of 
accuracy adequate for our purposes. In the following, we 
describe the quality criteria that we adopted to decide the 
inclusion or exclusion of a given galaxy from the analysis. 

Quality criterion W: weak sources 

The first criterion is based on the detection of a sufficient 
number of photons to obtain a reliable estimate of temper- 
ature and to perform the density deprojection. We require 
at least 100 counts in the region used for the temperature 
estimate and the presence of at least 5 bins in the sur- 
face brightness profile. For this reason, we discard 3C 075, 
3C 78, 3C 189, 3C 293, 3C 442, and UGC 7203. The counts 
in the X-ray image of UGC 7203, only marginally in excess 
of the background level, are insufficient to measure even 
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Table 2. Measured Bondi accretion rates and jet power. 
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Summary of the measured Bondi accretion rates and power of the FR I, CoreG, and Allen sample (marked with *). Column 
description: (1) name, (2) quality flag, (3) Bondi accretion radius in parsec, (4) ratio between the minimum radius reached by 
ClfANDRA resolution rmin and tb, (5) the logarithmic slope of the density profile, (6) temperature at the accretion radius 
[keV], (7) electron number density at the accretion radius [cm~"^], (8) Bondi accretion rate [MQyr~^], (9) predicted Bondi 
accretio n power [10^^ erg s~*], (10) radio jet power [10*^ erg s~*] (the error on Pj is Eissumed to be 0.4 dex, the dispersion of 
theE einz et aP l|2007t l relation linking it to the core power except for the Allen et al. objects for which we used their 
published values), (11) logarithm of the surface brightness at 1 kpc, Sikpc [counts s~^ arcsec"^]. 

its brightness profile. These 6 sources are marked with a Quality criterion C: complex morphologies 

W (weak) quality flag in Tabled 

In some cases, the X-ray galaxy morphology is quite com- 
plex, far away from the assumption of spherical symmetry 
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on which we based the deprojection technique. This is 
the case for three galaxies, namely 3C 084 (Perseus A), 
3C 317, and NGC 5128 (Cen A) as it can be seen from 
their respective Chandra images published in the litera- 
ture ( Fabian et al. I l2000t iBlantMT et al. raft et al 



These sources are marked with a C (Complex) qual- 
ity flag in Table El 

Quality criterion R: large rmin/i'B ratio 

The Bondi radius of most galaxies in our sample is typi- 
cally between 10 - 100 pc, which is always unresolved by 
Chandra data, with the only exception of UGC 7629 (aka 
NGC 4472). However, in some cases, the ratio between 
the minimum radius at which we are able to measure the 
electron density, Tmm, and tb is quite large, preventing 
a reasonable extrapolation of the density profile. For ex- 
ample, in 3C 348, the density at a Bondi radius of 34 pc 
should be obtained from a measurement at ~ 10 kpc. All 
galaxies (namely 3C 028, 3C 346, 3C 348, 3C 438, as well 
as the core galaxy UGC 9723) with rmin/i^B > 100 are 
excluded from the analysis at this level, and marked with 
a R flag. 

Quality criterion E: large extrapolation errors 

We required that the error in the extrapolation of the 
logarithm of the electron density to the accretion radius, 
as described in Sect. 13. 4[ must be smaller 1. We discarded 
for this reason 3C 31, 3C 264, and 3C 346 in the FR I 
sample and UGC 968, UGC 5902, UGC 7797, NGC 1316, 
and NGC 5419 among the CoreG (all marked with an E). 

We note that the errors in rig of the remaining sources 
are substantially smaller than the adopted threshold and 
its precise value is therefore not important. Similarly, we 
also verify that all sources excluded by the R criterion 
would also fall into this flag category, since the extrap- 
olation over a factor of 100 in radius always produces a 
large error in Ps , which indicate that the limit chosen for 
rmin/r_B IS not important. 

To summarize, we have 6 FR I sources and 8 CoreG for 
with an accurate estimate of the accretion power at the 
Bondi radius, reported in Table [21 to be added to the 9 
sources of the I Allen et al.l sample, for a total of 23 objects. 



4. The jet kinetic power 

Several estimators of jet kinetic power in radio galax- 
ies have been devised. A widely used method is based 
on the link between extended radio emission and jet 
power (jRawlings fc Saunderd Il99ll: IWillott et all Il999l ) 
that however suffers from substantial uncertainties (e.g. 
variations in environment, long-term time variability ef- 
fects), which are particularly severe when dealing with 
individual objects. Furthermore, the connection between 
extended radio emission and Pj is calibrated only for FR II 
radio galaxies and cannot be simply extrapolated to FR I. 



Results from X-ray studies with CHANDRA and 
XMM-Newton demonstrated the ability of a jet asso- 
ciated with an active galactic nucleus to blow cavities 
or 'bubbles' in the sur rounding X-ray emitting-gas (e.g 
Boehringer et al.lll993l : I Fabian et al.ll2000l : iForman et al 



20051 ). These surface brightness depressions associated 
with radio sources provide us with the most robust 
and precise estimati on of the jet kinetic power (e.g 
Churazov et al. I l2002h . assuming that the jet energy in- 
flates these bubbles at the sound speed. 

Measuring the pressure and density of the gas inside 
the X-ray bubbles, we can calculate the work done by 
the jet on the surrounding medium to inflate the cavi- 
ties {W (X pV). There are different sources of uncertainty 
in the measurements of cavity energy: the projected bub- 
ble size and the shape and composition of the contained 
plasma. The energy output is probably underestimated 
for adiabatic losses, cavity disruption, und etected cavi 
tics, and the omission o f shock energy (e.g. iNusser et al 
I2OO6; Bin nev et al.ll2007l ) . All these factors provide an un- 
certainty in the jet power of at least a factor of a few. 
Unfortunately, only a limited number of objects clearly 
shows surface brightness depressions associated with a ra- 
dio source. 



Heinz et al. 1 (I2OO7I) derived a relation between the ra- 
dio core luminosity L^, and the jet kinetic power Pj us- 
ing a sample of 15 galaxies that showed jet-driven X- 
ray cavities. The slope of the relation is determined by 
theoretical assumption about the jet emission, modeled 
by a power-law spectrum, while the normalization is ob- 
tained empirically. The possible drawbacks of this ap- 
proach are related to nuclear variability (as opposed to 
the time-averaged value obtained from the X-ray cavities) 
and Doppler boosting. Nonetheless, they showed that a 
signiflcant relation between these two quantities exists in 
the form: 



P. 



where Pq is a constant which has a value LOl^'g 10^^ 
erg s~^, when Lq is fixed at Lq = 7 10^^ ergHz^^ s^^. 
This relation then provides us with a simple and rather 
robust (with a scatter of 0.4 dex) estimate of the total 
jet power based only on the measure of the core power, a 
quantity available for all sources of our sample, with the 
only exception of 3C 028. 



We note that since the research bv I Heinz et al.l ([20071), 
two further measure ments of Pj based o n X-ray cavities 
have been published ( Shurkin et al.l[2007h . These two ob- 
jects (namely NGC 1399 and NGC 4649) closely follow 
the Pj - hcore relation, providing further support of its 
validity. Given the good agreement we do not need to re- 
estimate the value of Pq- 
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P. 


23 
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-17.6 
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Mbh 


P. 


44 


0.46 


0.48 


0.001 


2.05±0.20 


-17.4 


0.94 




P. 


43 


0.73 


0.72 


8x10"** 


1.17±0.11 
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0.67 



Col. (1) and (2) the correlated variables; Col. (3) the number of objects used; Col. (4) the linear correlation coefficient r; Col. 
(5) the Spearman's rank correlation p and (6) the probability of no correlation between the variables; Col. (7) the slope of the 
correlation; Col. (8) the intercept q; Col. (9) the rms scatte r from the correlat ion. * : Ps-Pj relation using only data from the 

9 objects in I Allen et al] l|2006l '). 




Fig. 4. Bondi accretion power versus jet power. The stars 
represent FR I objects, the circles ar e the CoreG and the 
triangles represent the objects from lAUen et al.l sample. 
We also plot the best-fit linear relation(solid line). 

5. Results 

5.1. The connection between jet and accretion power 

The estimated accretion rates at the Bondi radius in the 
galaxies of our sample are in the range from ~ 3 x 10^** — 
4 X 10"^ M0 year^^, corresponding to an available accre- 
tion power of ~ 2 X 10^^ — 2 x 10'*^erg s~^ 

In Fig. m we compare the estimated jet power and 
Bondi accretion power for the 23 objects for which we 
were able to measure this quantitjU. The Bondi accretion 
power and the jet power are strongly correlated, with a 



Spearman's rank coefficient of p = 0.76. The probability 
that this is produced by from a random distribution is 
only 3xl0~^. The linear correlation coefficient is instead 
r = 0.84 (see Table|3]for a summary of the statistical anal- 
ysis). The best- fit linear relations were derived to be the 
bisectrix of the fits using the two quantities as indepen 



dent v ariables; this followed the suggestion of llsobe et al 



With respect to lAllen et al.l we exchanged the variables on 
the two axis and used Pb ~ Mbc^ directly without any as- 
sumption about the efficiency of conversion of rest mass into 
energy. 



(119901 ) and is the preferable for problems that require the 
symmetrical treatment of the variables. The best-fit linear 
relation, weighted for the errors of both the quantities, is 

log(Pj) = (1.10±0.11)xlog(PB) - 1.91±0.20 

in 10^^ erg s^^ units. The rms dispersion of the data from 
the fit is 0.4 dex. Therefore, the jet power is, within the 
uncertainties, linearly related to the Bondi rate over 3 or- 
ders of magnitude with a scatter of less than a factor of 3. 
Unfortunately, all but 6 FR I had to be excluded from the 
analysis of the density profile (due partly to the presence of 
bright point sources and to a larger than average distance) 
and this prevents being able explore the behavior of the 
Pb vs Pj relation at the high end of radio luminosities 
of the sample. The efficiency of the conversion between 
accretion and j et power is 1.2lo;4 %. We note that the 
slope found by lAUen et all corresponds, in our notation, 
to 1.30±0.34 with a normalization of 2.2^0° % for Pj = 
erg s~^, consistent within 1 sigma confidence level 
with our values. 

Since the Bondi accretion rate depends on the black- 
hole mass, the central gas density and temperature, it is 
important to assess whether the relation between Pj and 
P B , discussed above, is driven mostly by only one of these 
quantities. To estimate these dependences, we compare 
the jet power with each of these variables separately. 

First, we found that there is no correlation between 
the jet power and the temperature at r^. Secondly, we 
considered how the electron density Ue at r^ affects the 
jet power. Fig. [5] (left panel) shows a weak correlation 
between Ug and the jet power with a Spearman's rank 
coefficient of p = 0.29 (with an associated probability of 
no correlation of 18 %) and a linear correlation coefficient 
of r = 0.34. The derived relation has a slope of 1.13±0.21. 
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Thirdly, we correlated Pj with the black-hole masses 
of the 23 objects for which we were able to estimate the 
accretion rate: Fig. [5] (middle panel) shows a correlation 
with a slope ~2 and a rms scatter of 0.66 dex. However, we 
can consider the measurements of mass for all 44 galaxies 
in the sample, regardless of the possibility to estimate Ps. 
In particular, this enables us to include the strongest ra- 
dio sources, radically changing the picture. The strength 
of the correlation between Mbh and Pj is strongly re- 
duced (see Table [3] and Fig. [5l right panel) using the en- 
tire sample. Most importantly, we find objects at a given 
black-hole mass with Pj differing by 3 orders of magni- 
tude and indeed the rms scatter of the residual is almost 
of 1 order of magnitude. T his result was presented in 
Capetti fc Balmaverdd ( 20051 ) where we discussed this is- 
sue using the same sample but without the further require- 
ment of the availability of Chandra data: while CoreG 
define a correlation between Mbh and Lcore (and conse- 
quently with Pj), FR I radio galaxies have much larger 
radio core luminosities at the similar value of M^^f which 
produces a large scatter in the Pj-MBH plane. 

We conclude that the jet power shows a connection to 
the Bondi accretion rate, which is stronger than those to 
Mbh, Ue, or T separately. This indicates that Pj vs Pb 
is the primary relation that does not descend from any of 
these quantities separately. 



5.2. Jet power and surface brightness profiles 

The analysis presented in the previous Section revealed 
the presence of a strong correlation between the jet power 
and the central properties of the galaxy, in particular with 
the gas density and black-hole mass. However, radio galax- 
ies of different power differ also ion terms of the global 
properties of their X-ray emitting hot coronae. 

These differences are clearly evident in a comparison 
of their surface brightness profiles (shown in Fig. [6] for 
the sub sample of 23 galaxies with estimated Pb)- Their 
SEP are described well by power-laws, with a relatively 
small rang e of slopes (<Alog5/Alogr> = -1.7±0.6j3 while 
they differ substantially in normalization. The remaining 
sources of the sample show a similar behavior with an even 
increased range of intensitj{f|. To quantify these differences, 
we estimated the value of the surface brightness at 1 kpc, 
(hereafter Sikpc) a radius accessible for almost all galaxies, 
by interpolating locally the profile with a power-law. The 
galaxies of our sample cover 3.5 orders of magnitude in 
^ikpc (see Tabled]). 

The intensity of the brightness profile at 1 kpc turns 
out to be strongly correlated with the radio-jet power (see 
Fig. El right panel). The correlation coefficient of this re- 
lation is 0.73, the rms is 0.67, and the Spearman rank co- 



efficient is 0.72. The probability, associated with this rank 
value, that this is produced by a random distribution is 
8x10"^. The slope of the correlation is 1.17±0.11. 

A significant advantage of this approach, based sim- 
ply on the surface brightness profile, is the possibility of 
including in the analysis also galaxies for which the data 
were of insufficient quality to estimate the electron den- 
sity at the Bondi radius. In particular, this enables us to 
extend the range of jet power probed to 4.5 orders of mag- 
nitudes. This is due to two reasons. As already noted, on 
the one hand, the most luminous radio galaxies did not, in 
general, have reliable density profiles due to the contam- 
ination of the nuclear source and to their average larger 
distance (e.g. 3C 84 and 3C 346). On the other hand, we 
find that the objects with the lowest radio luminosity cor- 
respond to hot coronae of lower surface brightness, which 
were then excluded from the analysis due to insufficient 
count rates in their X-ray images (e.g. UGC 5902 and 
UGC 9723). 

Therefore, the jet power of a given object can be pre- 
dicted with an rms uncertainty of only a factor of 5 by 
looking at the intensity of the X-ray emission produced 
by the galactic hot corona. 

Since the brightness profiles are directly related to the 
gas electron density, the large range of brightnesses ob- 
served reflects the different distributions of hot gas in the 
centers of these galaxies. A link between brightness profile 
and jet power is expected because the electron density is 
one of the critical parameters that, with the BH mass, de- 
termine the accretion flow. It must however be noted that, 
in the case of Bremsstrahlung emission, the expected de- 
pendence is in the form Pj oc iSifcpc"'^- U 

This suggests that our assumption of simple homol- 
ogy in the hot coronae, described only by the normal- 
ization of the brightness profile, is probably an over- 
simplification. Other variables can influence the relation 
between Pj and Sikpc, for example, a dependence of the 
SBP shape (e.g. its core radius or outer slope) on lumi- 
nosity. There is an indication of the existence of a fun- 
damental plane for the X-ray emitting gas, linking the 
temperature wi th the surface brightn ess and half-light (in 
X-ray) radius (jPiehl fc Statlej l2005l ). It is likely that a 
more robust parametrization of the hot corona properties, 
not simply its intensity at a fixed radius, might reveal an 
even stronger link with the AGN properties, but this is 
beyond the scope of this paper. 

An alternative explanation for the steeper than ex- 
pected dependence between brightness profi le and jet 



The average slope corresponds to the large scale behavior 
of an isothermal /3 model with /3 = 0.45 

^ Only the galaxy UGC 7203, showing extended emission 
only marginally in excess of the background level, is excluded 
from this analysis 



power originates in a result described by iBest et al 
( 20051 ). They found that the total mass of gas cooling 
within an isothermal hot corona is Mcooi oc Lx /T and 
that this has the same dependence on black-hole mass 
as the radio- loud AGN fraction Afcooi oc M^^. This im- 
plies that a linear relation exists between the cooling rate. 



® We tested that the inclusion of the effects on temperature 
of the count rates produces a marginal improvement in the 
correlation, but leaves the slope unchanged. 
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Fig. 5. The electron density at rs versus the jet power (left panel) for the sub sample of 23 objects for which we are 
able to measure the accretion power, with the linear best-fit linear relation. Middle panel) black-hole masses versus 
the jet power for this same sub sample. In the right panel, we compare the linear correlation (dashed line) obtained 
from the 23 objects (surrounded by a large circle) with that obtained from the entire sample studied (44 objects). 




Fig. 6. Left panel: surface brightness profiles of the 23 objects, zooming into the region around a radius of 1 kpc. 
They are described well by power laws, with a relatively small range of slopes, while they differ substantially in 
normalization. Right panel: Surface brightness at 1 kpc (Sikpc) in units of counts s^"'^ arcsec"^, versus the jet power 
for the extended sample of 43 galaxies. The solid line represents the best fit linear relation. 



radio-power, and the X-ray luminosity (and hence with 
the surface brightness) of the hot corona. 

Despite our inability to interpret in detail the corre- 
lation between Pj and Sikpc, the phenomenological con- 
nection between the brightness profile and jet power has 
several important ramifications. First of all, considering 
the galaxies with the highest radio luminosity, we already 
noted in Fig. O (right panel) that they show a substantial 
upward scatter in the plane comparing radio luminosity 
and black-hole mass. They effectively break the link be- 



tween hj and SMBH mass defined by the less luminous 
radio-sources. However, Fig. [HI (right panel) shows that 
these bright radio sources are associated with the bright- 
est coronae. Their higher radio luminosity, with respect to 
other galaxies of similar black-hole mass, is driven by the 
higher density of their hot ISM, producing a higher accre- 
tion rate and a more luminous extended X-ray emission. 

Furthermore, the results presented in Sect. 15.11 con- 
nect the jet power with the gas properties (and black-hole 
mass) at the Bondi radius. Conversely, Fig. [6] (right panel) 
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sets a link between the global properties of the hot coro- 
nae and the level of AGN activity. The implication of this 
result is that a substantial variation in the jet power in a 
given galaxy must be accompanied by a global change in 
its ISM properties. This can only occur when a dramatic 
event, such as a major merger, takes place. This is an in- 
dication that the accretion flow is stable over significant 
timescales and that, most likely, the long-term level of 
jet power remains stable over the interval between major 
mergers. 

Finally, the correlation between Pj and Sikpc provides 
us with a simple and rather robust method to estimate the 
jet power, with an accuracy of better than a factor of 5, 
based only on the measurement of the brightness of the hot 
gas at a fixed radius, 1 kpc for the sample considered here. 
Clearly, this is applicable to X-ray data of relatively poor 
quality and/or spatial resolution that would not enable 
a full spectral fitting and deprojection down to a spatial 
scale comparable to the Bondi radius. Indeed, this already 
found a useful application in the context of this work, as 
it enabled us to understand the properties of the bright 
PR I. We expect to be able to use this method also, for 
example, for more distant (and hence more powerful) radio 
galaxies. 

6. Constraints on the accretion process and on its 
radiative efficiency 

In Balmaverde fc Capettil ( 20061 ). we noticed that the 
X-ray nuclei of CoreG have luminosities of typically 
L/L^dd ^ 10^^ — 10^^ expressed in fraction of the 
Eddington luminosity. Purthermore, the tight correlations 
between radio, optical, and X-ray nuclear emission argue 
in favor of a jet origin for the nuclear emission. In PR I 
radio galaxies, the dominant contribution from jet emis- 
sion is confirm ed directly by the hi gh polarization of their 



process. This may be the case in the presence of an out- 
fiow (Ad yection Dominated Inflow /Ou tflow Solutions, or 
ADIOS, iBlandford fc Beerelmanl [l999 h or strong convec- 
tion (Convection Dominate d Accretion Plows, or CDAP, 
Quataert fc Gruzinov 2000l ) in which most gas circulates 
in convection eddies rather than accreting onto the black 
hole. These arguments apply also to the galaxies con- 
sidered here, since the available accretion powers deter- 
mined in the Sect. [3] exceed by several orders of magnitude 
the observed nuclear X-ray luminosities hx- More quan- 
tit atively, the comparison of the valu es of Ly reported 



in I Balmaverde fc Capetti ( 2006[ ) and IBalmaverde et al 
(200^ and Pb correspond to typical ratios of Lx / Ps 
- lO-'* - 10~^ 

In the case of the radio-loud objects considered here, 
we have further information that can be used to constrain 
the properties of the accretion process, related to the pres- 
ence of jets. In fact the mass flowing along the jet, mj, 
cannot exceed the mass rate at the radius of accretion 
disk, rhdisk- The jet power, for a matter-dominated jet, 
can be expressed as Pj ^ TrhjC^ where P is the jet Lorentz 
factor. If we denote with / the fraction of gas crossing 
the Bondi radius that actually reaches the accretion disk 
(that is TOdisk = /"t-b) and since the correlation between 
jet and Bondi power gives us 

Pj - 0.012 niBC^ 
the condition rrij < mdisk translates into a lower limit to 
/ as 

/ > 0.012 P-I 

Assuming the estimate of jet Lorentz for PR I ra dio galax- 
ies, P ~ 5, derived from lCiovannini et al.l (|2001l) this limit 
becomes 

/ > 0.0024 

Numerical simulations performed by IZanni et al. ( 2007 ) 
suggest that plasma within the accretion disk can be chan- 



optical nuclei (|Capetti et al.ll2007[ ). This imphes that the neled into the jet with a rather high efficiency. 



observed nuclear emission does not originate in the ac- 
cretion process and the values reported above should be 
considered as upper limits. 

These results add to the already vast literature report- 
ing emission corresponding to a very low fraction of the 
Eddington luminosity associated with accretion onto su- 
permassive black holes. This prompted the idea that in 
these objects accretion occurs not only at a low rate but 
also at a low radiative efficiency, suc h as in the Adyectio n 
Dominated Accretion Plows (ADAP. lNaravan fc YilllQQSl ). 
The ADAP models have been rather successful in model- 
ing the observed nuclear spectrum in severa l galaxies, such 
as e.g . the Galactic Cente r and NGC 4258 ( Naravan fc Yil 
19951 : iLasota et all Il996h . However, even ADAP mod- 



0.2 — 0.55 rhdisk from which / ^ 0.01. 

On the other hand, the correlation between the nuclear 
X-r ay luminosities and the radio -core luminosities derived 
in ( Balmaverde fc Capetti 20061 ) sets a limit to the disk 



Di Matteo et al.l l2000t 



els over-predict the observed em ission in the nuclei of 
nearby bright elliptical galaxies (ll 
Loewenstein et al. l l200lUPellegrinill2005h . 

This suggested the possibility that a substantial frac- 
tion of the mass included within the Bondi's accretion 
radius might not actually reach the central object, thus 
further reducing the radiative emission from the accretion 



emission as 

log{Lx) < lAO + log{iyLR) 

The X-ray luminosity can be expressed as 

Lx = VxtridiskC^ = VxfniBC^ = 0.012~^rixfPj 

where rj^ is the efficiency of conversion of rest mass into en- 
ergy radiated within the 2-10 keV band. Taking advantage 
of the link between radio-core luminosity and jet power 
discussed in Sect. IH we obtain: 

logLx = 24.84 + log{r]^f) + 12/17 logLn 

leading to 

logiVxf) < 5/17 logLR - 13.74 

The limit on rj^f is most stringent for the objects of lower 
radio luminosity. Por example, for Lr = lO"^"^ erg/s/Hz0 

^ corresponding to log uLr — 36.7 
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this corresponds to?73;/<1.6xl0 ^. When combined to 
the limit of / obtained above, this impUes 

Tj^ < 1.3 X 10""* r - 6.5 X 10""* 

To derive the total disk radiative efficiency we have to 
estimate the bolometric correction corresponding to the 2- 
10 keV band. This is a complex task, since little is known 
about the spectral energy distribution of our sources and, 
in addition, their broad -band nuclear s pectru m is dom- 
inated by jet emission. 'Ma rconi et"all (|2004l ) estimated 
bolometric corrections at different levels of luminosity by 
using suitably-built AGN templates and found that the 
2-10 keV bolometric correction decreases with decreas- 
ing luminosity, down to ~ 8 for a bolometric luminos- 
ity of ^ 10"*^ erg s~^. By adopting this value, we derive 
a bolometric disk radiative efficiency of 77 < 0.005. This 
limit is substantially smaller than the 'canonical' value of 
0.1 and of the estimates derived from matching the mass 
function of local SMBH with the AGN r elics, 77 - 0.08 
(jMarconi et all 12004 [Shankar et all 120071) . This is an in- 
dication that accretion in these objects occurs not only at 
a smaller rate, but also at lower efficiency than in standard 
accretion disks. 



7. Summary and conclusions 

We have presented our results concerning the connection 
between the accretion rate of hot gas, estimated in the 
Bondi's spherical approximation, and the jet power in a 
sample of low-luminosity radio galaxies. The sample was 
formed by combining radio galaxies with FR I morphology 
extracted from the 3C catalogue and early-type galaxies 
with a 'core' optical brightness profile. Galaxies of this lat- 
ter sub sample have been shown to represent "miniature" 
radio galaxies, being drawn from the same population of 
early-type galaxies as the FR I hosts and differing from 
"classical" FR I radio galaxies only in a lower level of nu- 
clear activity. This combined FR I/CoreG sample covers 
the full range of radio-luminosity of FR I, spanning five 
orders of magnitude. 

In the Chandra public archive, there were observations 
available for 44 objects of the sample. These data were 
used to derive profiles of X-ray brightness and of gas tem- 
perature that, once deprojected, provided a measurement 
of the electron gas density at the Bondi radius. When com- 
bined with the mass of the central black hole, these values 
provided us with an estimate of the accretion rate of hot 
gas, Pb- The quahty and/or spatial resolution of the data 
was sufficient to obtain a measure of Pb at a level of ac- 
curacy adequate for our purposes for a sub sample of 23 
galaxies. 

To estimate the jet power, we took advantage of an 
empirical relation linking the jet kinetic power, Pj, to the 
radio-core luminosity. This method freed us from the need 
of detecting X-ray cavities inflated by a radio jet to mea- 
sure its power and enabled us to estimate Pj with an rms 
error of 0.4 dex for all galaxies in our sample. 



We found that a tight, linear relation links the accre- 
tion with the jet power whose normalization sets the effi- 
ciency of the conversion between accretion and jet power 
to ~ 1.2%. These results stre ngthen and extend the valid- 
ity of the results obtained bv I Allen et al. I (|2006h based on 
the analysis of 9 low-luminosity radio galaxies. 

In addition, we found that the jet power is closely con- 
nected also with the surface brightness of the X-ray emis- 
sion, measured at the fixed radius of 1 kpc. This is a simple 
quantity to measure; it was therefore possible to consider 
in the analysis all galaxies, including those for which the 
data were of insufficient quality to estimate the Bondi ac- 
cretion rate. This method could be extended to objects of 
the highest radio luminosity, which did not provide use- 
ful estimates of Ps , showing that they are associated with 
the brightest coronae. This implies that the different levels 
of radio luminosity are not only associated with nuclear 
quantities (such as the gas density and black-hole mass) 
but also to global differences in the structure of the galac- 
tic hot coronae. A substantial change in the jet power must 
be accompanied by a global change in its ISM properties, 
driven e.g. by a major merger. 

At the estimated accretion rates at the Bondi radius, 
Mb , even ADAF models substantially over-predict the ob- 
served X-ray nuclear emission suggesting that a substan- 
tial fraction of the mass crossing the Bondi radius does 
not reach the central object. On the other hand, the mass 
rate at the outer radius of the accretion disk must be suffi- 
cient to maintain the mass flow along the jet. This requires 
that at least a mass fraction of > 0. 0024Mb reaches the 
accretion disk. The X-ray nuclear luminosity correspond- 
ing to this accretion rate does not exceed the observed 
luminosities of the nuclear X-ray sources (having assumed 
a bolometric correction of 8 for the 2-10 keV band) only 
when the radiative disk efficiency is 77 < 0.005, substan- 
tially smaller than the 'canonical' value of 0.1. This is an 
indication that accretion in these objects occurs not only 
at a smaller rate, but also at lower efficiency than in stan- 
dard accretion disks. 

The picture which emerges is that the accretion of hot 
gas associated with the host galaxy corona represents the 
dominant process of fueling for FR I radio galaxies across 
their full range of radio luminosity, spanning from the 
faintest level of detectable activity to radio-luminosities 
marking the transition between FR I and FR II objects. 
The differences of AGN power across this wide range of 
luminosity are determined by the available supply of gas, 
which we found to be closely related to the global prop- 
erties of the hot corona. These results provide a strong 
evidence for the presence of a feedback process linking the 
host galaxy with the active nucleus since the same hot gas 
that provides the supply to the supermassive black hole 
is also the repository of the kinetic energy liberated into 
the ISM by the jets. Clearly, it will be of great interest 
to explore how the hot gas content is connected to the 
evolution of the host galaxy and to its merger history. 

Another interesting issue is the fate of the gas within 
the Bondi accretion radius. While on the one handm the 



14 



Balmaverde et al.: The accretion mechanism in radio galaxies 



spherical accretion provides, a posteriori, a good descrip- 
tion of gas inflow in the central regions of these active 
galaxies, only a small fraction of the gas is allowed to 
reach the accretion disk. 

At small radii, the spherical inflow must therefore 
break and the bulk of the accreting gas must be removed 
from its large-scale motion toward the central black hole. 
The (almost) linear relation between jet power and accre- 
tion rate indicates that the efficiency of this mechanism 
must have only a very weak dependence on the amount 
of accreting gas. Apparently, also the fraction of gas ulti- 
mately launched into the relativistic jet is essentially in- 
dependent of the jet power. Accounting for this complex 
phenomenology represents an important challenge if we 
want to preceed in our understanding of the process of 
accretion and jet formation in active nuclei. 

A crucial question that must also be addressed is which 
manifestation of an active nucleus can be supported by hot 
accretion. Although the galaxies considered span a large 
range of radio power, it is important to establish whether 
hot accretion is sufficient to power FR II radio galaxies 
of even higher luminosity or, alternatively, above some 
threshold an additional source of fuel, possibly associated 
with an ISM phase of lower temperature. The results pre- 
sented here refer only to radio-loud AGN that channel a 
signiflcant fraction of the accretion power into jet kinetic 
energy. However, hot accretion might also be important in 
powering radio-quiet AGN. We will explore this possibility 
in a forthcoming paper. 

Acknowledgements. We acknowledge the usage of the 
Hyper Leda database (http:/ /leda.univ-lyonl.frj . 

References 

Allen, S. W., Dunn, R. J. H., Fabian, A. C, Taylor, G. B., & 

Reynolds, C. S. 2006, MNRAS, 372, 21 
Auriemma, C, Perola, G. C., Ekers, R. D., et al. 1977, A&A, 

57, 41 

Balmaverde, B. & Capetti, A. 2006, A&A, 447, 97 (BC06) 
Balmaverde, B., Capetti, A., & Grandi, P. 2006, A&A, 451, 35 
Bennett, A. S. 1962, MmRAS, 68, 163 
Best, P. N. 2004, MNRAS, 351, 70 

Best, P. N., Kauffmann, G., Heckman, T. M., et al. 2005, 

MNRAS, 362, 25 
Binney, J., Bibi, F. A., & Omma, H. 2007, MNRAS, 377, 142 
Birzan, L., Rafferty, D. A., McNamara, B. R., Wise, M. W., & 

Nulsen, P. E. J. 2004, ApJ, 607, 800 
Blandford, R. D. & Begelman, M. C. 1999, MNRAS, 303, LI 
Blanton, E. L., Sarazin, C. L., & McNamara, B. R. 2003, ApJ, 

585, 227 

Boehringer, H., Voges, W., Fabian, A. C., Edge, A. C., & 

Neumann, D. M. 1993, MNRAS, 264, L25 
Bondi, H. 1952, MNRAS, 112, 195 
Burns, J. O. 1990, AJ, 99, 14 

Capetti, A., Axon, D. J., Chiaberge, M., et al. 2007, A&A, 471, 
137 

Capetti, A. & Balmaverde, B. 2005, A&A, 440, 73 (CB05) 
Chiaberge, M., Capetti, A., & Celotti, A. 1999, A&A, 349, 77 
Churazov, E., Sunyaev, R., Forman, W., & Bohringer, H. 2002, 
MNRAS, 332, 729 



de Ruiter, H. R., Parma, P., Capetti, A., Fanti, R., & Morganti, 

R. 2002, A&A, 396, 857 
Di Matteo, T., Quataert, E., Allen, S. W., Narayan, R., & 

Fabian, A. C. 2000, MNRAS, 311, 507 
Diehl, S. & Statler, T. S. 2005, ApJ, 633, L21 
Donzelh, C. J., Chiaberge, M., Macchetto, F. D., et al. 2007, 

ApJ, 667, 780 

Fabian, A. C, Sanders, J. S., Ettori, S., et al. 2000, MNRAS, 
318, L65 

Forman, W., Nulsen, P., Heinz, S., et al. 2005, ApJ, 635, 894 
Giovannini, G., Cotton, W. D., Feretti, L., Lara, L., & Venturi, 

T. 2001, ApJ, 552, 508 
Gliozzi, M., Sambruna, R. M., & Brandt, W. N. 2003, A&A, 

408, 949 

Hardcastle, M. J., Evans, D. A., & Croston, J. H. 2007, 

MNRAS, 376, 1849 
Heinz, S., Merloni, A., & Schwab, J. 2007, ApJ, 658, L9 
Isobe, T., Feigelson, E. D., Akritas, M. G., & Babu, G. J. 1990, 

ApJ, 364, 104 

Kraft, R. P., Forman, W., Jones, C, et al. 2000, ApJ, 531, L9 
Kriss, G. A., Cioffi, D. F., & Canizares, C. R. 1983, ApJ, 272, 
439 

Lasota, J.-P., Abramowicz, M. A., Chen, X., et al. 1996, ApJ, 
462, 142 

Lauer, T. R., Ajhar, E. A., Byun, Y.-L, et al. 1995, AJ, 110, 
2622 

Loewenstein, M., Mushotzky, R. F., Angelini, L., Arnaud, 

K. A., & Quataert, E. 2001, ApJ, 555, L21 
Marconi, A. & Hunt, L. K. 2003, ApJ, 589, L21 
Marconi, A., Risaliti, G., GiUi, R., et al. 2004, MNRAS, 351, 

169 

Mauch, T. & Sadler, E. M. 2007, MNRAS, 375, 931 

Narayan, R. & Yi, 1. 1995, ApJ, 444, 231 

Nusser, A., Silk, J., & Babul, A. 2006, MNRAS, 373, 739 

Pellegrini, S. 2005, ApJ, 624, 155 

Quataert, E. & Gruzinov, A. 2000, ApJ, 539, 809 

Rawlings, S. & Saunders, R. 1991, Nature, 349, 138 

Shankar, F., Weinberg, D. H., & Miralda-Escude', J. 2007, 

ArXiv e-prints, 710 
Shurkin, K., Dunn, R. J. H., Gentile, G., Taylor, G. B., & 

Allen, S. W. 2007, ArXiv e-prints, 710 
Smith, E. P. & Heckman, T. M. 1989, ApJ, 341, 658 
Tremaine, S., Gebhardt, K., Bender, R., et al. 2002, ApJ, 574, 

740 

Valentijn, E. A. & Bijleveld, W. 1983, A&A, 125, 223 
Willott, C. J., Rawlings, S., BlundeU, K. M., & Lacy, M. 1999, 

MNRAS, 309, 1017 
Zanni, C, Ferrari, A., Rosner, R., Bodo, G., & Massaglia, S. 

2007, A&A, 469, 811 



Appendix A: 



Balmaverde et al.: The accretion mechanism in radio galaxies 



15 







X 




i 




I 








z 

z 




I 

i 










X 




i w 

i 




I 








I 




i 

i 






' i ■ 






X 




I I 




I I 




















I X 




X 

X ^ 










1 ^ i 


X 




I 


I 




a 


I 







16 



Balmaverde et al.: The accretion mechanism in radio galaxies 



I 




: s 




i i 




I 




X 








i 








i 




i 




X 

1 




i 








I 

ii J 












I 









ill 



T ' ^ 

I i i 






Balmaverde et al.: The accretion mechanism in radio galaxies 




Fig. A.l. The X-ray surface brightness profiles of the 43 galaxies of the sample, having excluded UGC 7203. 



18 



Balmaverde et aL: The accretion mechanism in radio galaxies 




3.10 1.00 



100.00 1000.00 0.01 



Radius [Kpc] 

3C 189 



Radius [Kpc] 
3C 296 



Radius [Kpc] 

3C 438 



10.00 0.01 



100.00 0.01 



0.10 1.00 



Radius [Kpc] 

3C 078 



Radius [Kpc] 
3C 264 



Rddius [Kpc] 

3C346 



Radius [Kpc] 
3C 449 



■10,00 0.01 



Radius [Kpc] 
3C 83.1 



Radius [Kpc] 
3C 270 



Rodijs [Kpc] 

3C 348 



100.00 0.01 0.10 



100.00 1000.00 



1.000 10.000 100.000 



Balmaverde et al.: The accretion mechanism in radio galaxies 



19 



Rodus [Kpc] 

UGC 5297 



'00.00 0.01 



UGC 7386 



- ^ 0,1000 r 



G.OOO 0.01 



Radius [Kpc] 

UGC 8745 



Radij. [Kpc] 

NGC 1316 



10. OD 0.01 



Radius [Kpc] 

NGC 3557 



-E ^ 0.0100 r 



10.00 0,01 



10,00 0,01 



Fig. A. 2. The electron density profiles of the 30 objects, whose data quality is sufficient to perform the deprojection 
of the brightness profile (see Sect. 13. 5p . The solid lines represent the best-fit power-law relation to the density profile. 



